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Abstract 

The structure of low-energy collective states in proton-dehcient N = 28 isotones is analyzed 
using structure models based on the relativistic energy density functional DD-PC1. The relativistic 
Hartree-Bogoliubov model for triaxial nuclei is used to calculate binding energy maps in the f3 - 
7 plane. The evolution of neutron and proton single-particle levels with quadrupole deformation, 
and the occurrence of gaps around the Fermi surface, provide a simple microscopic interpretation 
of the onset of deformation and shape coexistence. Starting from self-consistent constrained energy 
surfaces calculated with the functional DD-PC1, a collective Hamiltonian for quadrupole vibrations 
and rotations is employed in the analysis of excitation spectra and transition rates of 46 Ar, 44 S, and 
42 Si. The results are compared to available data, and previous studies based either on the mean- 
field approach or large-scale shell- model calculations. The present study is particularly focused on 
44 S, for which data have recently been reported that indicate pronounced shape coexistence. 

PACS numbers: 21.10.-k, 21.60.Jz, 21.60.Ev 



I. INTRODUCTION 



Shapes of neutron-rich nuclei far from stability have extensively been explored in many 
experimental and theoretical studies. The evolution of ground-state shapes in an isotopic 
or isotonic chain, for instance, is governed by changes of the underlying shell structure 
of single-nucleon orbitals. In particular far from the /^-stability line, the energy spacings 
between single-nucleon levels change considerably with the number of neutrons or protons. 
This can lead to reduced spherical shell gaps, and in some cases spherical magic numbers 
may partly or entirely disappear The reduction of spherical shell closure often leads to 
the occurrence of ground-states deformation and, in a number of cases, to the coexistence 
of different shapes in a single nucleus. 

In recent years a number of studies have been devoted to the investigation of the fragility 
of the N = 28 magic number in neutron- rich nuclei j^J. In /3-stable nuclei the Z or N = 28 
shell closure is the first magic number produced by the spin-orbit part of the single-nucleon 
potential, which lowers the fj/2 orbital with respect to the P3/2 and thus forms a spherical 
shell gap at nucleon number 28. However, as a number of experimental investigations have 



shown 3Ml3l|. in the proton-deficient iV = 28 isotones below 48 Ca the spherical shell gap is 
progressively reduced and the low-energy spectra of 46 Ar, 44 S, and 42 Si display evidence of 
ground-state deformation and shape-coexistence. 

Both large-scale shell model (SM) calculations (4-6, 9 
(SCMF) models 



2, B SI 



19] and self-consistent mean-field 



26] have been employed in the theoretical description of these 
phenomena. The basic advantages of the SM approach include the ability to simultaneously 
describe all spectroscopic properties of low-lying states, the use of effective interactions that 
can be related to microscopic inter-nucleon forces, and the description of collective properties 
in the laboratory frame. On the other hand, since SM effective interactions depend on the 
choice of active shells and truncation schemes, there is no universal shell-model interaction 
that can be used for all nuclei. 

A variety of structure phenomena, including regions of exotic nuclei far from the line of /In- 
stability and close to the nucleon drip-lines, have been successfully described with mean-field 
models based on the Gogny interaction, the Skyrme energy functional, and the relativistic 



meson-exchange effective Lagrangian 



The SCMF approach to nuclear structure 



enables a description of the nuclear many-body problem in terms of a universal energy den- 
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sity functional (EDF). When extended to also take into account collective correlations, this 
framework provides a detailed microscopic description of structure phenomena associated 
with shell evolution. Compared to the SM, the strong points of the mean-field approach are 
the use of global functionals, the treatment of arbitrarily heavy systems, model spaces that 
include all occupied states (no distinction between core and valence nucleons, no need for 
effective charges) and the intuitive picture of intrinsic shapes. 

A quantitative description of shell evolution, and in particular the treatment of shape 
coexistence phenomena, necessitates the inclusion of many-body correlations beyond the 
mean-field approximation. The starting point is usually a constrained Hartree-Fock plus 
BCS (HFBCS), or Hartree-Fock-Bogoliubov (HFB) calculation of the binding energy surface 
with the mass quadrupole components as constrained quantities. In most studies calculations 
have been restricted to axially symmetric, parity conserving configurations. The erosion of 
spherical shell-closures in nuclei far from stability leads to deformed intrinsic states and, in 
some cases, mean-field potential energy surfaces with almost degenerate prolate and oblate 
minima. In order to describe nuclei with soft potential energy surfaces and/or small energy 
differences between coexisting minima, it is necessary to explicitly consider correlation effects 
beyond the mean-field level. The rotational energy correction, i.e. the energy gained by the 
restoration of rotational symmetry, is proportional to the quadrupole deformation of the 
intrinsic state and can reach several MeV for a well deformed configuration. Fluctuations of 
quadrupole deformation also contribute to the correlation energy. Both types of correlations 
can be included simultaneously by mixing angular momentum projected states corresponding 
to different quadrupole moments. The most effective approach for configuration mixing 
calculations is the generator coordinate method (GCM), with multipole moments used as 
coordinates that generate the intrinsic wave functions. 

In recent years several accurate and efficient models, based on microscopic energy density 
functionals, have been developed that perform restoration of symmetries broken by the static 
nuclear mean field, and take into account quadrupole fluctuations. However, while GCM 
configuration mixing of axially symmetric states has routinely been employed in structure 
studies, the application of this method to triaxial shapes presents a much more involved 
and technically difficult problem. Only the most recent advances in parallel computing 
and modeling have enabled the implementation of microscopic models, based on triaxial 
symmetry-breaking intrinsic states that are projected on particle number and angular mo- 
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mentum, and finally mixed by the generator coordinate method 3CH33]. 

In an approximation to the full GCM for five-dimensional quadrupole dynamics, a col- 
lective Hamiltonian can be formulated that restores rotational symmetry and accounts for 
fluctuations around mean-field minima. The dynamics of the five-dimensional Hamilto- 
nian for quadrupole vibrational and rotational degrees of freedom is governed by the seven 
functions of the intrinsic deformations /3 and 7 : the collective potential, the three vibra- 
tional mass parameters, and three moments of inertia for rotations around the principal 
axes. These functions are determined by microscopic mean-field calculations using a uni- 
versal nuclear EDF. Starting from self-consistent single-nucleon orbitals, the corresponding 
occupation probabilities and energies at each point on the constrained energy surfaces, the 
mass parameters and the moments of inertia are calculated as functions of the deformations 
j3 and 7. The diagonalization of the resulting Hamiltonian yields excitation energies and 



collective wave functions t 



rat can be used to calculate various observables, such as electro- 



magnetic transition rates 34j, |35| . In this work we employ a recent implementation of the 
collective Hamiltonian for quadrupole degrees of freedom in a study of shape coexistence 
and low-energy collective states in N = 28 isotones. 

Both non-relativistic and relativistic energy density functionals have been used in SCMF 
studies of the erosion of the N = 28 spherical shell gap. One of the advantages of using 
relativistic functionals, particularly evident in the example of N = 28 isotones, is the natural 
inclusion of the nucleon spin degree of freedom, and the resulting nuclear spin-orbit potential 
which emerges automatically with the empirical strength in a covariant formulation. In 



the present analysis we use the new relativistic functional DD-PC1 [36]. Starting from 
microscopic nucleon self-energies in nuclear matter, and empirical global properties of the 
nuclear matter equation of state, the coupling parameters of DD-PC1 were fine-tuned to the 
experimental masses of a set of 64 deformed nuclei in the mass regions A pa 150 — 180 and 
A pa 230 — 250. The functional has been further tested in calculations of medium-heavy and 
heavy nuclei, including binding energies, charge radii, deformation parameters, neutron skin 
thickness, and excitation energies of giant multipole resonances. The present calculation of 
N = 28 isotones, therefore, presents an extrapolation of DD-PC1 to a region of nuclei very 
different from the mass regions where the parameters of the functional were adjusted, and 
thus a test of the global applicability of DD-PC1. 

Section [TT] includes a short review of the theoretical framework: the relativistic Hartee- 
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Bogoliubov model for triaxial nuclei, and the corresponding collective Hamiltonian for 
quadrupole degrees of freedom. The evolution of shapes in the N = 28 isotones is ana- 
lyzed in Sec. IHIt the quadrupole constrained energy surfaces determined by DD-PC1, and 
the resulting low-energy collective spectra, in comparison to available data and previous 
SCMF and SM calculations. Section IIVI summarizes the results and ends with an outlook 
for future studies. 



II. THEORETICAL FRAMEWORK 



A. 3D relativistic Hartee-Bogoliubov model with a separable pairing interaction 



The relativistic Hartee-Bogoliubov model [28|, |29j provides a unified description of 
particle-hole (ph) and particle-particle (pp) correlations on a mean-field level by combining 
two average potentials: the self-consistent mean field that encloses long range ph correla- 
tions, and a pairing field A which sums up pp-correlations. In the present analysis the 



in the 
is 



mean-field potential is determined by the relativistic density functional DD-PC1 36 
ph channel, and a new separable pairing interaction, recently introduced in Refs. 37 
used in the pp channel. 

In the RHB framework the mean-field state is described by a generalized Slater deter- 
minant |$) that represents the vacuum with respect to independent quasiparticles. The 
quasiparticle operators are defined by the unitary Bogoliubov transformation, and the cor- 
responding Hartree-Bogoliubov wave functions U and V are determined by the solution of 
the RHB equation. In coordinate representation: 



hp — m — A A 

-A* -h* D + m + X 





(1) 



In the relativistic case the self-consistent mean-field corresponds to the single-nucleon Dirac 
Hamiltonian hp, m is the nucleon mass, and the chemical potential A is determined by the 
particle number subsidiary condition such that the expectation value of the particle number 
operator in the ground state equals the number of nucleons. The pairing field A reads 

A ab (r, r') = - V ab cd(r, r')«od(r, r'). (2) 

c,d 
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where V a bcd{v, r') are the matrix elements of the two-body pairing interaction, and the indices 
a, b, c and d denote the quantum numbers that specify the Dirac indices of the spinor. 
The column vectors denote the quasiparticle wave functions, and E k are the quasiparticle 
energies. 

The single-particle density and the pairing tensor, constructed from the quasiparticle 
wave functions 

Pcd (r,r') = J2VcUr)V dk (r'), (3) 

fc>0 

K cd (r,r>) = J2u: k (r)V dk (r'), (4) 

fc>0 

are calculated in the no-sea approximation (denoted by k > 0): the summation runs over 
all quasiparticle states k with positive quasiparticle energies E k > 0, but omits states that 
originate from the Dirac sea. The latter are characterized by quasiparticle energies larger 
than the Dirac gap (» 1200 MeV). 

In most applications of the RHB model the pairing part of the Gogny force j39] was 
used in the particle-particle (pp) channel. A basic advantage of the Gogny force is the 
finite range, which automatically guarantees a proper cut-off in momentum space. However, 
the resulting pairing field is non-local and the solution of the corresponding Dirac-Hartree- 
Bogoliubov integro-differential equations can be time-consuming, especially for nuclei with 
non-axial shapes. For that reason a separable form of the pairing interaction was recently 



. The interaction 



introduced for RHB calculations in spherical and deformed nuclei [3j 
is separable in momentum space: (k\V So \k'} = -Gp(k)p(k') and, by assuming a simple 
Gaussian ansatz p(k) = e _a2fc2 , the two parameters G and a were adjusted to reproduce 
the density dependence of the gap at the Fermi surface in nuclear matter, calculated with 



a Gogny force. For the D1S parameterization of the Gogny force 39], the corresponding 
parameters of the separable pairing interaction take the following values: G = —728 MeVfm 3 
and a = 0.644 fm. When transformed from momentum to coordinate space, the force takes 
the form: 

V(r t , r 2 , ri, r' 2 ) = G5(R- R') P{r)P{r') 1 - (1 - P ff ) , (5) 

where R = \{ri + r 2 ) and r = r\ — r 2 denote the center-of-mass and the relative coordi- 
nates, and P(r) is the Fourier transform of p(k): 



(Ana 2 ) 



The pairing interaction is of finite range and, because of the presence of the factor 5 (R — R), 
it preserves translational invariance. Even though 5 (R — R') implies that this force is not 
completely separable in coordinate space, the corresponding pp matrix elements can be rep- 
resented as a sum of a finite number of separable terms in the basis of a three-dimensional 
(3D) harmonic oscillator. The interaction of Eq. (jSJ) reproduces pairing properties of spher- 
ical and axially deformed nuclei calculated with the original Gogny force, but with the 
important advantage that the computational cost is greatly reduced. 

To describe nuclei with general quadrupole shapes, the Dirac-Hartree-Bogoliubov equa- 
tions ((!]) are solved by expanding the nucleon spinors in the basis of a 3D harmonic oscillator 
in Cartesian coordinates. In the present calculation of N = 28 isotones complete conver- 
gence is obtained with NJ 1 ^ = 10 major oscillator shells. The map of the energy surface as 
a function of the quadrupole deformation is obtained by imposing constraints on the axial 
and triaxial quadrupole moments. The method of quadratic constraint uses an unrestricted 
variation of the function 

(#>+£ Cv((Q 2/1 )-g 2M ) 2 , (7) 

11=0,2 

where (H) is the total energy, and {Qiy) denotes the expectation value of the mass 
quadrupole operators: 

Q 20 = 2z 2 -x 2 -y 2 and Q 22 = x 2 - y 2 . (8) 
is the constrained value of the multipole moment, and C 2m the corresponding stiffness 



constant 



401]. 



B. Collective Hamiltonian in Five Dimensions 

The self-consistent solutions of the constrained triaxial RHB equations, i.e. the single- 
quasiparticle energies and wave functions for the entire energy surface as functions of the 
quadrupole deformation, provide the microscopic input for the parameters of a collective 



Hamiltonian for vibrational and rotational degrees of freedom [3J]. The five quadrupole 
collective coordinates are parameterized in terms of the two deformation parameters (3 and 
7, and three Euler angles (<ft, 6, ip) = Q, which define the orientation of the intrinsic 
principal axes in the laboratory frame. 

-f^coll = 2~Ub + Trot + Koll > (9) 



with the vibrational kinetic energy: 



T v ih — — 



h 2 



+ 



1 



d_ 

dp 



(3 sin 37 
and rotational kinetic energy: 



d_ 



v d d i v d 

r did pF d 

-^3 7j B /37 - + --^-sin3 7 ^- 



rot 



T 2 

k=i K 



(10) 



(11) 



V^ou is the collective potential. J k denotes the components of the angular momentum in 
the body- fixed frame of a nucleus, and the mass parameters Bpp, Bp^, _B 77 , as well as the 
moments of inertia Ik, depend on the quadrupole deformation variables j3 and 7: 



l k = 4B k (3 2 sin 2 (7 - 2kn/3) 



(12) 



Two additional quantities that appear in the expression for the vibrational energy: r = 
B1B2B3, and w = BppB^ — B 2 ^, determine the volume element in the collective space. 

The dynamics of the collective Hamiltonian is governed by the seven functions of the 
intrinsic deformations (3 and 7: the collective potential, the three mass parameters: Bpp, 
Bp*,, .B 77 , and the three moments of inertia Ik- These functions are determined by the 
microscopic nuclear energy density functional and the effective interaction in the pp channel. 
The moments of inertia are calculated from the Inglis-Belyaev formula: 



Ik 



t=l,2,8, 



Ei + Ej 



(13) 



where k denotes the axis of rotation, the summation runs over proton and neutron quasiparti- 
cle states \ij) = |$), and |<3>) represents the quasiparticle vacuum. The mass parameters 
associated with the two quadrupole collective coordinates q = (Q20) and q 2 = (Q22) are 
calculated in the cranking approximation: 

h 2 



l (i) 



J [IV 



where 



M (n) „Mo, q 2 ) = L (E~+E~y> 



(14) 



(15) 
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Finally, the potential V co u in the collective Hamiltonian Eq. (Q is obtained by subtracting the 
zero-point energy (ZPE) corrections from the total energy that corresponds to the solution 
of constrained RHB equations, at each point on the triaxial deformation plane 34]. 

The Hamiltonian Eq. (jH]) describes quadrupole vibrations, rotations, and the coupling of 
these collective modes. The corresponding eigenvalue problem is solved using an expansion 
of eigenfunctions in terms of a complete set of basis functions that depend on the deforma- 



tion variables /3 and 7, and the Euler angles 0, 9 and i/j 



34|. The diagonalization of the 



Hamiltonian yields the excitation energies and collective wave functions: 

K M i^i^) = E V£*G0,7)*i«c(n)- (is) 

KeAI 

The angular part corresponds to linear combinations of Wigner functions 



*mk(P) = J ig^* 1 ^) + (-I)^m-k(^)] , (17) 

and the summation in Eq. ffTBT) is over the allowed set of the K values: 

f 0, 2, .... I for I mod 2 = 
AI = < (18) 
I 2,4, 1 for / mod 2 = 1 . 

Using the collective wave functions Eq. (fT6l) . various observables can be calculated and 
compared with experimental results. For instance, the quadrupole E2 reduced transition 
probability: 

5(E2; al a'l') = ^^|<aT||^(i<;2)||a/}| 2 , (19) 

where M.(E2) is the electric quadrupole operator, local in the collective deformation vari- 
ables. 



III. EVOLUTION OF SHAPES IN THE N=28 ISOTONES 



A. Quadrupole binding energy maps 

The 3D relativistic Hartree-Bogoliubov model, with the functional DD-PC1 in the 
particle-hole channel and a separable pairing force in the particle-particle channel, enables 
very efficient constrained self-consistent triaxial calculations of binding energy maps as func- 
tions of quadrupole deformation in the (3 — 7 plane. The resulting single-quasiparticle en- 
ergies and wave functions provide the microscopic input for the GCM configuration mixing 

9 



of angular-momentum projected triaxial wave functions, or can be used to determine the 
parameters of the collective Hamiltonian for vibrations and rotations: the mass parameters, 
the moments of inertia, and the collective potential. The solution of the corresponding 
eigenvalue problem yields the excitation spectra and collective wave functions that are used 
in the calculation of electromagnetic transition probabilities. This approach is here applied 
to the low-energy quadrupole spectra of N = 28 isotones. 
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FIG. 1: (Color online) Self-consistent RHB triaxial quadrupole constrained energy surfaces of 
N = 28 isotones in the (3 — 7 plane (0 < 7 < 60°). For each nucleus energies are normalized with 
respect to the binding energy of the global minimum. The contours join points on the surface with 
the same energy (in MeV). 

Figure [1] displays the self-consistent RHB triaxial quadrupole constrained energy surfaces 
of N = 28 isotones in the /3 — 7 plane (0 < 7 < 60°), calculated using the DD-PC1 energy 
density functional, plus the separable pairing force Eq. (JSJ) in the particle-particle channel. 
For each nucleus energies are normalized with respect to the binding energy of the absolute 
minimum. The contours join points on the surface with the same energy. 

Starting from the spherical doubly-magic 48 Ca, we consider the even-even N = 28 isotones 
obtained by successive removals of proton pairs. The binding energy maps display a rich 
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variety of rapidly evolving shapes, and clearly demonstrate the fragility of the N = 28 shell. 
By removing a pair of protons from 48 Ca, the energy surface of the corresponding isotone 
46 Ar becomes soft both in (3 and 7, with a shallow extended minimum along the oblate 
axis. Only four protons away from the doubly magic 48 Ca, DD-PC1 predicts a coexistence 
of prolate and oblate minima at (/?, 7) = (0.34,0°) and (0.27,60°), respectively, in 44 S. The 
two minima are separated by a rather low barrier of less than 1 MeV and, therefore, one 
expects to find pronounced mixing of prolate and oblate configurations in the low-energy 
collective states of this nucleus. For 42 Si the binding energy displays a deep oblate minimum 
at (/3,7) = (0.35,60°), whereas a secondary, prolate minimum is calculated ~ 2.5 MeV 
higher. Finally, with another proton pair removed, the very neutron-rich nucleus 40 Mg 
shows a deep prolate minimum at (/?, 7) = (0.45,0°). 



We note that similar binding energy surfaces were also obtained in recent studies |41l. 142 1 
based on the self-consistent Hartree-Fock-Bogoliubov (HFB) model, using the finite-range 
and density- dependent Gogny D1S interaction. On the mean-field level the only qualitative 
difference is found for 40 Mg. For this nucleus the present calculation predicts a saddle 
point on the oblate axis, whereas a secondary local oblate minimum is obtained in the HFB 
calculation with the Gogny force. 

The variation of mean-field shapes in an isotopic, or isotonic, chain is governed by the 
evolution of the underlying shell structure of single-nucleon orbitals. The formation of 
deformed minima, in particular, can be related to the occurrence of gaps or regions of low 
single-particle level density around the Fermi surface. In Figs. [2] - [5] we plot the neutron 
and proton single-particle energy levels in the canonical basis for 46 Ar, 44 S, 42 Si, and 40 Mg, 
respectively. Solid (black) curves correspond to levels with positive parity, and (red) dashed 
curves denote levels with negative parity. The dot-dashed (blue) curves correspond to the 
Fermi levels. The neutron and proton levels are plotted as functions of the deformation 
parameters along closed paths in the (3 — 7 plane. The panels on the left and right display 
prolate (7 = 0°) and oblate (7 = 60°) axially-symmetric single-particle levels, respectively. 
In the middle panel of each figure the neutron and proton levels are plotted as functions of 
7, for a fixed value of the axial deformation \/3\ at the approximate position of the mean-field 
minima: \0\ = 0.2 for 46 Ar, \p\ = 0.3 for 44 S, and \p\ = 0.4 for 42 Si and 40 Mg. In this way, 
starting from the spherical configuration, we follow the single-nucleon levels on a path along 
the prolate axis up to the approximate position of the minimum (left panel), then for this 
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fixed value of \(3\ the path from 7 = 0° to 7 = 60° (middle panel) and, finally, back to the 
spherical configuration along the oblate axis (right panel). Negative values of beta denote 
axial deformations with 7 = 60°, that is, points along the oblate axis. 

Neutron 




0.0 0.1 20 40 60 -0.1 0.0 



p Y (deg) p 

Proton 




P Y (deg) p 

FIG. 2: (Color online) Single-neutron and single-proton energy levels of 46 Ar, as functions of the 
deformation parameters along closed paths in the f3 — 7 plane. Solid (black) curves correspond 
to levels with positive parity and (red) dashed curves denote levels with negative parity. The 
dot-dashed (blue) curves corresponds to the Fermi levels. The panels on the left and right display 
prolate (7 = 0°) and oblate (7 = 60°) axially-symmetric single-particle levels, respectively. In the 
middle panel of each figure the neutron and proton levels are plotted as functions of 7, for a fixed 
value of the axial deformation |/3| at the approximate position of the mean-field minimum. 

Figures [2] - [5] elucidate the principal characteristics of structural changes in neutron-rich 
iV = 28 nuclei: the near degeneracy of the and S1/2 proton orbitals, and the reduction 
of the size of the N = 28 shell gap [2j. Between the doubly magic 48 Ca and 46 Ar the 
spherical gap A^ = 28 decreases from 4.73 MeV to 4.48 MeV (cf. Table [I]), in excellent 
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FIG. 3: (Color online) Same as described in the caption to Fig. [2] but for the nucleus 44 S. 

agreement with data: from 4.80 MeV in 48 Ca to 4.47 MeV in 46 Ar Nevertheless, the 

gap between occupied and unoccupied neutron levels in 46 Ar is still largest at the spherical 
configuration, as shown in the upper panel of Fig.[2j We note, in particular, the agreement of 
the calculated energies of spherical neutron states with experimental single-neutron energies 
For the proton states shown in the lower panel, the largest gap is found at \/3\ = 0.2 and 
7 = 60°, that is, on the oblate axis. The competition between the spherical configuration 
favored by neutron states and the oblate shape favored by proton states, leads to the shallow 
extended oblate minimum shown in Fig. [TJ Two protons less, and the spherical N = 28 gap 
is reduced by another 620 keV to 3.86 MeV in 44 S. The largest gap between neutron states 
is not the spherical one like in 46 Ar, however, but at the oblate deformation \/3\ ~ 0.3 and 
7 = 60° (upper panel of Fig. [3]). The removal of two protons lowers the energy of the 
corresponding Fermi level, and for 44 S the largest gap is found on the prolate axis (lower 
panel of Fig. |3]). The formation of the oblate neutron and prolate proton gaps is at the 
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FIG. 4: (Color online) Same as described in the caption to Fig. [2] but for the nucleus 42 Si. 



origin of the coexistence of deformed shapes in 44 S (cf. Fig. [T]). In 42 Si both neutron and 
proton gaps are on the oblate axis resulting in the pronounced oblate minimum at \/3\ ~ 
0.35. Finally, the deep prolate minimum at (3 ~ 0.35 in 40 Mg arises because of the neutron 
gap and, especially pronounced, proton gap on the prolate axis. We note that the largest 
neutron gap for this nucleus is still on the oblate side but, because the protons strongly 
favor the prolate configuration, it produces only a saddle point on the oblate shown 
in Fig. [TJ 

The erosion of the spherical N = 28 shell is also shown in Table HI where we include the 
DD-PC1 RHB theoretical neutron N = 28 spherical energy gaps, and the corresponding 
values of the axial deformation for the minima of the quadrupole binding energy maps 
of 48 Ca, 46 Ar, 44 S, 42 Si, and 40 Mg. Both experiment and theory point toward a strong 
reduction of the N = 28 gap as more protons are removed and, thus, the isotones become 
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FIG. 5: (Color online) Same as described in the caption to Fig. [2] but for the nucleus 40 Mg. 



more neutron-rich. N=28 is the first "magic" number produced by the spin-orbit part of 
the single-nucleon potential and, therefore, a relativistic mean-field model automatically 
reproduces the N=28 gap because it naturally includes the spin-orbit interaction and the 
correct isospin dependence of this term, as it was already shown in the axial RHB calculation 



of neutron-rich N=28 nuclei 



22j. Experimentally, indirect evidence of the erosion of the gap 



has been obtained by following the evolution of excitation energies of the 2f state and the 




E2 transitions in iV = 28 isotones and neighboring nuclei [3H5|, |9|, [13(. The experimental 
results can be reproduced by both mean-field 22, 24] and shell model [isl calculations. As 
shown in Table [H the DD-PC1 RHB calculation predicts a reduction of the spherical iV = 28 
shell gap from 4.73 MeV in the doubly-magic nucleus 48 Ca to 2.03 MeV in the well-deformed 
40 Mg. We note that the theoretical values of the spherical shell gap for 48 Ca and 46 Ar are 
very close to data: 4.80 MeV in 49 Ca, and 4.47 MeV in 47 Ar, obtained by neutron stripping 
reactions jsl, 
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TABLE I: The DD-PC1 RHB theoretical neutron N = 28 spherical energy gaps, and the corre- 
sponding values of the axial deformation for the minima of the quadrupole binding energy maps 
of 48 Ca, 46 Ar, 44 S, 42 Si, and 40 Mg. Negative values of (3 denote oblate shapes. 
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FIG. 6: (Color online) The spectrum of 46 Ar calculated with the DD-PC1 relativistic density 



functional (left), compared to data 44] (right) for the excitation energy of 2f, and the reduced 
electric quadrupole transition B(E2) (in units of e 2 fm 4 ). The prediction for the electric monopole 
transition strength p 2 (i?0;0^ — > Of) is also included in the theoretical spectrum. 

B. Low-energy collective spectra 

Starting from constrained self-consistent solutions of the RHB equations, that is, using 
single-quasiparticle energies and wave functions that correspond to each point on the en- 
ergy surfaces shown in Fig. [lj the parameters that determine the collective Hamiltonian: 
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FIG. 7: (Color online) Same as described in the caption to Fig. [6] but for the nucleus 44 S. The 
data are from Refs. P. [b]]. 
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data are from Ref. 



11] 



Same as described in the caption to Fig. [6] but for the nucleus 42 Si. The 



the mass parameters B@p, B^, -B 77 , three moments of inertia as well as the zero-point 
energy corrections, are calculated as functions of the quadrupole deformations (3 and 7. The 
diagonalization of the resulting Hamiltonian yields the excitation energies and reduced tran- 
sition probabilities. In Figs. [6]- El we display the spectra of 46 Ar, 44 S, and 42 Si calculated 
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with the DD-PC1 relativistic density functional plus the separable pairing force Eq. (jSJ), in 
comparison to available data for the excitation energies, reduced electric quadrupole tran- 
sition probabilities B(E2) (in units of e 2 fm 4 ), and the electric monopole transition strength 
p 2 (E0; Of — > Of). We emphasize that this calculation is completely parameter-free, that is, 
by using the self-consistent solutions of the RHB single-nucleon equations, physical observ- 
ables, such as transition probabilities and spectroscopic quadrupole moments, are calculated 
in the full configuration space and there is no need for effective charges. Using the bare value 
of the proton charge in the electric quadrupole operator, the transition probabilities between 
eigenstates of the collective Hamiltonian can directly be compared to data. 
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FIG. 9: (Color online) Evolution of the characteristic observables E{2f) and B(E2; 2f — > Of) (in 
e 2 fm 4 ) with proton number in N = 28 isotones. The ratio between the excitation energies of the 
first 4 + and 2 + states is also displayed in the inset. The microscopic values calculated with the 
energy density functional DD-PC1 are shown in comparison with available data. 



Before considering the excitation spectra of individual nuclei and, in particular, shape 
coexistence in 44 S, in Fig. |H]we illustrate the evolution with proton number of characteristic 
collective observables: the excitation energy of the first 2 + state, the ratio E(4f)/E(2f), 
and the reduced transition probability B(E2; 2+ — > Of). The rapid decrease of the ratio 
E(4f)/E(2f) from w 2.8 in 40 Mg to ~ 2.1 in 48 Ca is characteristic for a transition from 
a deformed rotational nucleus to a spherical vibrator. Note, however, that even in the 
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case of 40 Mg the value of E(4f)/E(2f) is considerably below the rigid-rotor limit of 3.3. 
The excitation energy of the first excited state 2f can directly be compared to data. The 
calculated E(2f) increases with proton number toward the doubly magic 48 Ca, but the 
predicted rise in energy is not as sharp as in experiment. In fact, one expects that in 
deformed nuclei, e.g 42 Si, the calculated E(2f) is above the experimental excitation energy, 
because of the well-known fact that the Inglis-Belyaev formula Eq. ({TBI) predicts effective 
moments of inertia that are smaller than empirical values. The moments of inertia can 
generally be improved by including the Thouless-Valatin (TV) dynamical rearrangement 



contributions 



4l|, but the calculation of the TV moments of inertia 45j has not yet been 



implemented in the collective Hamiltonian used in the present calculation. The panel on the 
right of Fig. IH] displays the evolution with proton number of another characteristic collective 
observable: B(E2;2f — > Of) (in e 2 fm 4 ). The calculation reproduces the empirical decrease 
of B(E2; 2f — > Of) with proton number and, in particular, we notice the excellent agreement 
between the parameter-free theoretical predictions and data for 44 S and 46 Ar. 

Figure [6] displays the low-energy spectrum of 46 Ar. The excitation energy E[2f) is cal- 
culated considerably above the experimental state, whereas the B(E2; 2f — > Of) reproduces 
the experimental value. In the present analysis we particularly focus on 44 S, for which data 
that indicate shape coexistence were reported recently 13j. Already the data on the low 
energy of the first 2 + state and the enhanced B(E2;2f — > Of) of 63(18) e 2 fm 4 [5] pointed 
towards a possible deformation of the ground state of 44 S. More recently, the structure of 
this nucleus was studied by using delayed 7 and electron spectroscopy, and new data were 
reported for the reduced transition probability B(E2;2f — > Of) = 8.4(26) e 2 fm 4 , and the 
monopole strength p 2 (E0;0f ->• 0^) = 8.7(7) x 10" 3 Q- From a comparison to shell 
model calculations, a prolate-spherical shape coexistence was inferred, and a two-level mix- 
ing model was used to extract a weak mixing between the two configurations. The spectrum 
of 44 S calculated in this work is compared to available data in Fig. [71 The model nicely re- 
produces both the excitation energy and the reduced transition probability B[E2\ 2f — > Of) 
for the first excited state 2f, and the theoretical value for B(E2; Of — > 2f) is also in good 
agreement with data. The experimental ratio B(E2\ 2f 0f)/B(E2; 2f — > Of) is 7.5, and 
the calculated value is 5.2. The excitation energy of the state 0^~, however, is calculated 
much higher than the experimental counterpart. Together with the fact that the calculated 
monopole transition strength p 2 (E0; Of — > 0^)(xl0 3 ) = 23 is larger than the corresponding 
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experimental value of 8.7(7), this result indicates that there is more mixing between the 
theoretical states Of and than what can be inferred from the data. 

The low- lying state with the excitation energy 1.365 MeV, the rather weak inter- 
band transition probability B(E2; 2f — > Of) =8.4(26) e 2 fm 4 , and the monopole strength 
p 2 (E0; Of — > Of) = 8.7(7) x 10 -3 have been regarded as fingerprints of shape coexistence in 



44 c 



131 ]. One reason for the more pronounced mixing between the calculated 0^ and Of in 



this work and, consequently, the higher excitation energy of 0^~, could be the particular choice 
of the energy density functional and/or the treatment of pairing correlations 46) . The pre- 
dicted barrier between the prolate and oblate minima (cf. Fig. [T]) could, in fact, be too low. 
Another reason for the high excitation energy of 0^" could be the approximation used in the 
calculation of mass parameters (vibrational inertial functions). In the current version of the 
model the mass parameters are determined by using the cranking approximation Eqs. ( JT4"]) 
and (Tl5|) . in which the time-odd components (the so-called Thouless-Valatin dynamical re- 
arrangement contributions) are omitted. Recently an efficient microscopic derivation of the 
five-dimensional quadrupole collective Hamiltonian has been developed, based on the adia- 
batic self-consistent collective coordinate method 1471. In this model the moments of inertia 



and mass parameters are determined from local normal modes built on constrained Hartree- 
Fock-Bogoliubov states, and the TV dynamical rearrangement contributions are treated 
self-consistently. For the illustrative case of 68 Se, it has been shown that the self-consistent 
inclusion of the time-odd components of the mean-field can lead to an increase of the mass 
parameters by 30% ~ 200%, depending on the deformation. In fact, in the present calcula- 
tion an enhancement of the cranking masses by a factor ~ 2 brings the calculated excitation 
energies, and also the monopole strength p 2 (E0; Of — > Of), in very close agreement with the 
experimental spectrum. 

In Table [TT] we compare the experimental excitation energies of the states 2f, Of, and 
2f, the reduced transition probabilities B(E2; 2f ->• Of) (e 2 fm 4 ), B(E2; 2f ->• Of), and the 
monopole strength p 2 (E0;0f — > Of) x 10 3 in 44 S, to the results of the present work, the 
five-dimensional GCM(GOA) calculation with the Gogny D1S interaction J42], the angular- 
momentum projected GCM calculation restricted to axial shapes (AMP GCM) with the 
Gogny D1S interaction [25] , and to shell-model calculations 13] . One might notice that 
all three models based on constrained self-consistent mean-field calculations of the binding 
energy maps (curves in the case of axially-symmetric AMP GCM), reproduce the data with 
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similar accuracy. It is interesting that only the axially-symmetric calculation reproduces 
the very low excitation energy of the state 0^, whereas the result of the five- dimensional 
GCM(GOA) calculation, although it was also based on the Gogny D1S interaction, is even 
above the energy obtained with DD-PC1. Table HT1 shows that the best overall agreement 
with data is obtained in the shell- model (SM) calculation of Ref. [3], using the effective 



interaction SDPF-U 



181 ] for Oftw SM calculations in the sd — pf valence space, and with a 



particular choice of the proton and neutron effective charges. 



TABLE II: Excitation energies (in MeV) of the states 2^, Oj, and 2j, B{E2;2\ -> 0^) (e 2 fm 4 ) 
5(£2;2t" -> Oj), and the monopole strength p 2 (E0;0% -> 0^) x 10 3 in 44 S. The experimental 
values [y, [p| are compared to the results of the present work, the five-dimensional GCM(GOA) 
calculation with the Gogny D1S interaction 42|, the angular-momentum projected GCM calcula- 
tion restricted to axial shapes (AMPGCM) with the Gogny D1S interaction [251 ] , and to shell-model 
calculations fl^ ]. 
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mi) 




1.329(1) 


1.491 


1.267 


1.410 


1.172 


E(0t) 




1.365(1) 


2.852 


3.611 


1.070 


1.137 


mi) 




2.335(39) 


2.851 


2.557 


1.830 


2.140 


B(E2;2i -» 




63(18) 


72 


105 


75 


75 


B(E2;2+ -» 


o 2 + ) 


8.4(2.6) 


14 


6.3 




19 




)(xl0 3 ) 


8.7(7) 


23 


5.4 







Based on the data included in Table [TT] and on the SM calculation with the SDPF- 
U effective interaction, in Ref. 13J] it was deduced that 44 S exhibits a shape coexistence 



between a prolate ground state (/3 ~ 0.25) and a rather spherical state. The sequence of 
ground-state band states Of, 2+, 4£> an d ®i> * s connected by strong E2 transitions, and the 
excited states are characterized by the intrinsic quadrupole moment Qq ~ 60 e fm 2 . This 
sequence was interpreted as a rotational band of an axially deformed prolate shape with 
f3 ~ 0.25. The calculated 2i state has a smaller quadrupole moment Qo = —0.3 e fm 2 , 
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compatible with a spherical shape, and is connected by a strong E2 transition to the 



13|- 



state. These SM results, therefore, indicate a prolate-spherical shape coexistence in 44 S 

To analyze configuration mixing in the low-energy spectrum based on the functional DD- 
PC1, in Fig. [TUlwe plot the probability density distributions for the three lowest states of 
the ground-state band: 0^, 2f, and 4f, the state 0^", and the two states 2\ and 2jj~. For 
a given collective state Eq. ( TT6|) . the probability density distribution in the (/?, 7) plane is 
defined by: 

Pia{P,i)= E l^(A7)| 2 /3 3 , (20) 
KeAi 

with the normalization: 

r-2n 



I fid? Pi a {fin) |sin3 7 | d 7 = l. (21) 
Jo 

The probability distribution of the ground state 0+ indicates a deformation |/3| > 0.3, ex- 
tended in the 7 direction from the prolate configuration at 7 = to the oblate configuration 
at 7 = 60°. The average deformation is ((/?), (7)) = (0.32,26°), and the 7-softness reflects 
the ground-state mixing of configurations based on the prolate and oblate minima of the 
potential (cf. Fig. [T|). With the increase of angular momentum in the ground-state band, 
e.g. 2f, 4f, etc., the states are progressively concentrated on the prolate axis. For instance, 
((/3), (7)) = (0.35,23°) for 2±. The average /3-deformation in the ground-state band grad- 
ually increases because of centrifugal stretching. Again we note that the empirical value 
B(E2; 2^ — > O2 ) is accurately reproduced by the present calculation using just the bare 
proton charge. In contrast to the SM prediction 13] . the state 0^~ is predominantly prolate, 
although one notices a relatively large overlap between the wave functions of the states 0i 
and 0^". The mixing between these states is probably one of the reasons for the high exci- 
tation energy of the second + state, as predicted by the present calculation (cf. Fig. [7j). 
The probability distribution of the state 2jj~ is concentrated on the prolate axis, and this 
state is connected by a strong transition to 0^: B{E2\2\ — > 0^) = 66 e 2 fm 4 , comparable 
to B{E2\2\ — > 0^). We note, however, that for the "coexisting" band based on 0^" the 
calculated ratio E(A + )/E{2 + ) is only 2.33. 

The calculated second 2 + state displays a probability distribution extended in the 7- 
direction and peaked on the oblate axis. As shown in Fig. [7] and Table UH this state is 
very close to the experimental candidate for the 2% state, which was suggested to be at 



2335(39) keV by placing the 988 keV transition [6] on top of the 0^" or 2\ state 13J. The 
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FIG. 10: (Color online) Probability distribution Eq. (|20p in the j3 — 7 plane for the lowest collective 
states of 44 S, predicted by DD-PC1 energy density functional. 

theoretical 2^~ state can be interpreted as the (quasi)-7 band-head according to the strong 
E2 transitions to the states 3± and 4 J. For the three lowest 2 + states, in Table IIHI we 
include the percentage of the K = and K = 2 components in the corresponding collective 
wave functions Eq. ( TT6l) (K denotes the projection of the angular momentum on the intrinsic 
3-axis), as well as the spectroscopic quadrupole moments. The wave functions of the states 
2f and 2g" are dominated by K = components, and the spectroscopic quadrupole moments 
are negative (prolate configurations) with comparable magnitudes. The positive quadrupole 
moment of 2^ points to a predominant oblate configuration, and the ~ 80% contribution of 
the K = 2 component in the wave function confirms that this state is the band-head of a 
(quasi) 7-band (note the formation of the doublet 3+ and 4^). 

Finally, Fig. [8] shows the low-energy collective spectrum of 42 Si. Even though the excita- 
tion spectrum and transition pattern appear to be similar to that of 44 S, with the exception 
of a considerably weaker E2 transition 0^ — > 1\ (cf. Fig. [7]), the ground-state band of 
this nucleus is in fact based on the oblate minimum shown in the binding energy map of 
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TABLE III: Percentage of the K = and K = 2 components (projection of the angular momentum 
on the body-fixed symmetry axis) for the collective wave functions of the three lowest 2 + states in 
44 S, and the corresponding spectroscopic quadrupole moments (in e fm 2 ). 
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FIG. 11: (Color online) Same as described in the caption to Fig. [10] but for the nucleus 42 Si. 

Fig. [TJ This is nicely illustrated in Fig. [TT] where, just like in the case of 44 S in Fig. [TUJ, 
we plot the probability distributions of the collective wave functions 0^, 2^, and 4^, the 
state 0^~, and the two states l\ and 2jj~. The wave functions of the yrast states 0^, 2^, and 
4+ are concentrated along the oblate axis. The state is strongly prolate deformed, with 
a peak in the probability distribution at (3 ~ 0.5. This state has a much smaller overlap 
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with 2^ than in the case of 44 S, and this explains the correspondingly weaker transition. 
For 42 Si, therefore, the solution of the collective Hamiltonian based on the DD-PC1 func- 
tional, predicts a coexistence of the oblate yrast band and the prolate sequence built on the 
strongly deformed state Oj. As already shown in Fig. [HI the present calculation does not 



reproduce the exceptionally low excitation energy of the state 2\: 770(19) keV [ill ]. It is 



interesting, however, that the calculated excitation energy of this state is very close to the 



SM prediction obtained using the SDPF-NR effective interaction [TJJ. Only by removing 
from the SDPF-NR a schematic pairing Hamiltonian in the pf shell, that is, by using the 
new effective interaction SDPF-U jl8|. the 2 + excitation energies of the silicon isotopes can 
be brought in agreement with experiment. 

IV. SUMMARY 

Structure phenomena related to the evolution of single-nucleon levels and shells in 
neutron-rich nuclei present a very active area of experimental and theoretical research. 
Among the microscopic models that can be used for a theoretical analysis of these phe- 
nomena, the framework of nuclear energy density functionals (EDFs) presently provides 
a complete and accurate description of ground-state properties and collective excitations 
across the entire chart of nuclides. In this work we have used the recently introduced rela- 
tivists EDF DD-PC1 36] to study the erosion of the N = 28 spherical shell in neutron-rich 
nuclei and the related phenomenon of shape evolution and shape coexistence in the N = 28 
isotones 46 Ar, 44 S, 42 Si, and 40 Mg. Pairing correlations have been taken into account by em- 
ploying an interaction that is separable in momentum space, and is completely determined 
by two parameters adjusted to reproduce the empirical bell-shaped pairing gap in symmetric 



nuclear matter 38]. 

The N=28 shell closure is the first neutron magic number produced by the spin-orbit part 
of the single-nucleon potential and, therefore, a relativistic mean-field model automatically 
reproduces the N=28 spherical gap because it naturally includes the spin-orbit interaction 
and the correct isospin dependence of this term, as it was shown more than ten years ago 
in the axial RHB calculation of neutron- rich N=28 nuclei 22J. In particular, in the RMF 



approach there is no need for a tensor interaction to reproduce the isospin dependence 
(quenching) of the spherical N=28 gap in neutron-rich nuclei, as also shown in the present 
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work in Table [U compared to available data. 

The functional DD-PC1 was adjusted exclusively to the experimental masses of a set of 
64 deformed nuclei in the mass regions A w 150 — 180 and A pa 230 — 250. The present 
study of the N = 28 isotones thus presents an extrapolation of DD-PC1 to a completely 
different region of the nuclide chart, and a further test of the universality of nuclear EDFs. 
It is not at all obvious that such an extrapolation will produce results in agreement with 
experiment, especially in a detailed comparison with spectroscopic data. The fact that it 
does is remarkable, and justifies the approach to nuclear structure based on universal energy 
density functionals. 

Starting from self-consistent binding energy maps in the (3 — 7 plane, calculated in the 
relativistic Hartree-Bogoliubov (RHB) model based on the functional DD-PC1, a recent 
implementation of the collective Hamiltonian for quadrupole vibrations and rotations has 
been used to calculate the excitation spectra and transition rates of 46 Ar, 44 S, 42 Si, and 
40 Mg. The parameters that determine the collective Hamiltonian: the vibrational inertial 
functions, the moments of inertia, and the zero-point energy corrections, are calculated 
using the single-quasiparticle energy and wave functions that correspond to each point on 
the self-consistent RHB binding energy surface of a given nucleus. The diagonalization of 
the collective Hamiltonian yields the excitation energies and wave functions used to calculate 
various observables. 

The calculation performed in this work has shown that the relativistic functional DD- 
PC1 provides an accurate microscopic interpretation of the strong reduction of the iV = 28 
spherical energy gap in neutron-rich nuclei, and a quantitative description of the evolution of 
shapes in N = 28 isotones in terms of single-nucleon orbitals as functions of the quadrupole 
deformation parameters and 7. In particular, the predicted values for the spherical shell 
gap in 48 Ca (4.73 MeV) and in 46 Ar (4.48 MeV), are very close to the data: 4.80 MeV 
in 49 Ca and 4.47 MeV in 47 Ar. The solutions of the collective Hamiltonian based on DD- 
PC1 reproduce the evolution with proton number of characteristic collective observables the 
excitation energy of the first 2 + state, the ratio £'(4 ] f )/ E(2~l), and the reduced transition 
probability B(E2; 2+ — y Of). In the present work we have focused on 44 S, for which recent 
data point towards a coexistence of shapes with different deformations in the low-energy 
excitation spectrum. It has been shown that the formation of the oblate neutron and 
prolate proton gaps, illustrated in Fig. [3J is at the origin of the predicted coexistence of 
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deformed shapes in 44 S. The spectroscopic results have been compared to available data, to 
triaxial (collective Hamiltonian) and axial (generator coordinate method) calculations based 
on the Gogny D1S HFB self-consistent mean-field energy maps, and to recent shell-model 
(SM) calculations using the new SDPF-U effective interaction. The present results are in 
qualitative agreement with previous calculations based on the Gogny D1S HFB model and, 
in particular, reproduce the data on both the excitation energy of the first excited state 
2\ and the reduced transition probability B(E2;2f — > Of), and the theoretical value for 
B[E2; Of — > 2f) is also in good agreement with data. The experimental ratio B(E2; 2f — > 
0f)/B(E2;2f — > Of) is 7.5, and the calculated value is 5.2. The theoretical monopole 
transition strength p 2 (E0; Of — > 0^)(xl0 3 ) = 23 is somewhat larger than the corresponding 
experimental value of 8.7(7). The calculation of transition rates in the collective Hamiltonian 
model is completely parameter- free. One might notice that the results predicted by the 
functional DD-PC1 have been compared to those obtained using effective interactions that 
were fine-tuned to data that include also this mass region, or adjusted exclusively to data 
in this region of the mass table (shell- model interactions). The fact that a global density 
functional can even compete in a spectroscopic calculation with shell-model interactions 
specifically customized to this mass region, and the level of agreement with experiment, 
presents a valuable result. 

A discrepancy with respect to experiment in 44 S is the high excitation energy predicted 
for the state 0^~, a factor of two compared to data. It appears that the model predicts too 
much mixing between the two lowest + states, and this leads to an enhancement of the 
corresponding monopole transition strength. The pronounced mixing between the calculated 
0^ and states, and the resulting repulsion, could be at the origin of the high excitation 
energy of 0^. The most obvious reason is that this is an intrinsic prediction of the functional 
DD-PC1. To check this one would have to perform calculations using different functionals 



481 ] . However, since also the Gogny D1S + 5DCH model yields a similar result, the functional 



itself probably is not the main problem. A more probable reason is that the mass parameters 
calculated in the cranking approximation are simply too small, as discussed in Sec. IIIIBl 
Finally, the excited + state could also have pronounced non-collective components that 
are not included in our model space (2-quasiparticle contributions). This is certainly a 
possibility, and it would partially explain why the calculated B(E2) to the first 2 + state is 
larger than the experimental value. The shell-model calculation of Ref. [3] predicts the 
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excitation energy of 0^ in better agreement with experiment, but the calculated B(E2) for 
the transition from the first 2 + state is more than a factor two larger than the experimental 
value (only about 50% larger in the present calculation). Therefore, it appears that the 
structure of the second + state in 44 S remains an open problem. 

This present analysis of low-energy spectra of N = 28 isotones has clearly demonstrated 
the advantages of using EDFs in the description of deformed nuclei: an intuitive mean-field 
interpretation in terms of coexisting intrinsic shapes and the evolution of single-particle 
states, spectroscopic calculations performed in the full model space of occupied states, and 
the universality of EDFs that enables their applications to nuclei in different mass regions, 
including short-lived systems far from stability. 
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